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', Abstract 

o 

' An ideal gas of twodimensional Dirac fermions in the background of a pointlike 

■ magnetic vortex with arbitrary flux is considered. We find that this system acquires 

fractional electric charge at finite temperatures and determine the functional depen- 
dence of the thermal average and quadratic fluctuation of the charge on the tempera- 
ture, the vortex flux, and the continuous parameter of the boundary condition at the 



Oh! 
(D . 

r~| , location of the vortex. 

> : 

^ : I Introduction 

Spontaneous breakdown of continuous symmetries can give rise to topological defects (texture 
solitons) with rather interesting properties. A topological defect in threedimensional space, 
which is characterized by the nontrivial second homotopy group, is known as a magnetic 
monopole [H E], see also genuine Ref.jB]. Vacuum fluctuations of quantized Dirac fields 
result in the monopole becoming a CP symmetry violating dyon, i.e. acquiring nonzero 
(and fractional) electric charge [H El EI- More recently the effect of thermal fluctuations 
of quantized Dirac fields in the presence of the monopole has been considered, yielding the 
temperature dependence of the induced charge 13 El E] • 

A topological defect in twodimensional space, which is characterized by the nontrivial 
first homotopy group, is a cross section of the Abrikosov-Nielsen-Olesen magnetic vortex 
Pm ITT] . The vortex defect is described in terms of a spin-0 field which condenses and a 
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spin-1 field corresponding to the spontaneously broken gauge group; the former is coupled 
to the latter in the minimal way with constant econd- Single-valuedness of the condensate 
field and finiteness of the vortex energy imply that the vortex flux is related to econd'- 

$ = ^ /rfxV(x) = ^, (1.1) 

where V(x) is the vector potential of the spin-1 field, and the integral is over a path enclosing 
once the vortex tube. The quantized fermion field is coupled minimally to the spin-1 field 
with constant e - the elementary charge; thus, quantum effects depend on the value of 
e$. The case of Ccond = 2e (e$ = 1/2) is realized in ordinary Bardeen-Cooper-Schrieffer 
superconductors where the Cooper pair field condenses and, in addition, there are normal 
electron (pair-breaking) excitations. It remains still to be elucidated, whether other values 
of e$ are realized in nature, although there are claims that vortices with fractional e$ 7^ 1/2 
exist in chiral superfiuids and chiral and two-gap superconductors [121 HH] ■ 

The aim of the present paper is to consider the effect of thermal fluctuations of quantized 
Dirac fields ^ in the presence of the vortex defect with arbitrary value of e$, which results 
in the vortex acquiring fractional electric charge; the zero-temperature effect was considered 
earlier ^Tl^JUHI- Since continuous symmetry is not spontaneously broken at the core of 
the defect, it seems reasonable to exclude the region of the defect and to impose a boundary 
condition for quantized fields at the edge of this region. Thus, quantum effects depend both 
on e$ and real continuous quantity which parameterizes the most general varieties of 
boundary conditions (for more details see next Section). This setup should not be confused 
with the setup when fermions are quantized in the presence of an extensive magnetic field 
with finite flux and the region of the nonvanishing fleld strength is not excluded. The 
induced charge in the latter case was considered in Refs. [201 |2I1 1221 (zero temperature) and 
Refs. |23| I25j (nonzero temperature), and we shall compare the results of both setups in 
Section V. 

The operator of the second-quantized fermion fleld in a static background can be pre- 
sented in the form 



xi/(x, t) = ^e-^^'{^\E, X)aEx + ^ e-'^*(x|i?, A)6+, , (1.2) 

E>0 E<0 

where a'^^^ and {b~^x ^ex) are the fermion (antifermion) creation and destruction 
operators satisfying anticommutation relations, 

[(^EX, at;,y] ^ = [bEX, &^,A'] + = , A|^', A') , (1.3) 

and {x.\E, A) is the solution to the stationary Dirac equation, 

H{x\E,X) = E{^\E,X) , (1.4) 
H is the Dirac Hamiltonian, E is the energy and A is the set of other parameters (quan- 

■■■■■■.„, ■.„...■.■..,, „■., t „ - ... 

^This may be relevant for various particle physics models with applications ranging from early Universe 
cosmology to hot nuclear matter phenomenology, and even for condensed matter models, because effectively 
quasirelativistic fermions arise, in particular, in d-wave type II superconductors (see, e.g., Refs. [T^ 1151 ) 
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integration (with a certain measure) over continuous values of all quantum numbers. Con- 
ventionally, the operators of dynamical invariants are defined as bilinears of the fermion field 
operators, and, thus, comprizing: 

the energy operator (temporal component of the energy-momentum vector), 

po=iy d''x(^[^-',dt^]_-[dt^+,^]_) , (1.5) 

and the fermion number operator, 

N=^jd'^x [^+,^]_ , (1.6) 

where d is the space dimension. Operators ()1.5|) and ()1.6|) commute and are thus diagonal 
in the fermion and antifermion creation and destruction operators. 

The thermal average of the fermion number operator over the canonical ensemble is 
defined as (see, e.g., Ref. [26]) 

Sp exp(-/3P0) 

where T is the equilibrium temperature, is the Boltzmann constant, and 5*^ is the trace or 
the sum over the expectation values in the Fock state basis created by operators in Eq. (|1.3|) . 
Appropriately, the electric charge of the quantum fermionic system in thermal equilibrium 
is given by expression 



oo 

Q(T) = e{N) = -\ j dET{E) tanh Q/?^) 



-oo 



where the last equality is obtained by transforming the right hand side of Eq. ()1.7j) into an 
integral over the spectrum of the Dirac Hamiltonian (see, e.g., Ref. |23j). and the spectral 
density of the Dirac Hamiltonian (or density of states) is 

where Tr is the trace of an integro-differential operator in functional space: Trll = 
= J d'^xtr{x.\U\x.); tr denotes the trace over spinor indices only; note that the functional 
trace should be regularized and renormalized by subtraction, if necessary. 

Similarly, one gets expression for the quadratic fluctuation of the electric charge: 



oo 



f! / dE ^^^L. (1.10) 

4 J cosh' {lf3E) ^ ' 



Evidently, if the quadratic fluctuation becomes nonvanishing, then the corresponding dy- 
namical invariant ceases to be a sharp quantum observable. 

In the present paper we shall find electric charge ()1.8|) and its fluctuation p.lO|) in the 
d = 2 quantum fermionic system in the background of a single static topological defect which 
is a twodimensional cross section of the magnetic vortex. 
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II Self-adjointness of the Dirac Hamiltonian in the back- 
ground of the pointhke vortex defect 

The Dirac Hamiltonian in external magnetic field takes form 

i7 = - ieV(x)] + 7°m . (2.1) 

In 2 + 1-dimensional space-time (x, t) = {x^, x'^,t), the Clifford algebra has two inequivalent 
irreducible representations which can be differed in the following way: 

i^o^i^2 s = ±l. (2.2) 

Choosing the 7*^ matrix in the diagonal form, 

7° = ^^3 , (2.3) 

one gets 

where Ui, a2, and are the Pauli matrices, and Xi ^-^id X-i ^-^e the parameters varying in 
interval < Xs < 27r to go over to the equivalent representation. Note also that in odd- 
dimensional space-time the m parameter in Eq. ()2.H) can take both positive and negative 
values; a change of sign of m corresponds to going over to the inequivalent representation. 
A solution to stationary Dirac equation ()1.4j) with Hamiltonian ()2.H) can be presented 

as 

/ f (r FA pi-^V+iXs \ 



where polar coordinates = y (2^^)^ + (a;^)^ and ip = arctan(x^/x^) are introduced, and Z 
is the set of integer numbers. The magnetic field strength and its flux in the units of 27r are 
given by expressions (compare with Eq. ljl.H) ): 

5(x) = 9xV(x), $ = — /rf^x5(x). (2.6) 

27r J 

Since a single defect is considered, a support of the magnetic field strength is localized 
in a certain, let it be central, region of twodimensional space. It is evident that different 
functions -B(x) can give the same value of In general, a solution to the Dirac equation in 
an external magnetic field depends on the field configuration in two ways: there is a direct, 
or local, impact of -B(x) on the solution at the same point x (similar to the action of the 
classical Lorentz force), and there is an indirect (through vector potential) influence of the 
field strength on the behaviour of the solution in regions out of the field strength's support 
(similar to the quantum-mechanical Bohm-Aharonov effect j27|). Namely the latter effects 
are those which interest us, since, as it has been already noted in Introduction, the central 
region (i.e. the region of the field strength's support) is excluded, and at its edge a boundary 
condition is imposed on the solution to the Dirac equation. 

It might be anticipated that the solution in the outer region depends on the flux rather 
than the local features of the field strength in the central region. However, a closer look at 
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the situation when the central region is not excluded suggests that the local features of the 
field strength may influence the behaviour of the solution out of the central region. Really, 
changes of a profile of the field strength influence strongly the behaviour of the solution in 
the central region, and this, due to the continuity and smoothness properties of the solution 
as a solution to a differential equation, entails changes in its behaviour in the outer region. 
Thus, when the central region is excluded, our purpose is not to stick to a limited set of 
boundary conditions but, instead, to extend this set maximally in order to cover all possible 
types of the behaviour of the solution near the boundary and, perhaps, all plausible profiles 
of the field strength in the excluded region. 

How to achieve this purpose in general, remains to be a question. However, a recipe is 
available under a simphfying assumption that finite size of the excluded region is neglected: 
then the most general conditions are those ensuring self-adjointness of the Dirac Hamiltonian, 
and they are labelled by self-adjoint extension parameter (see, e.g., Ref.|2Hl)- Although 
is physically interpreted as the CP violating vacuum angle in the d = 3 case of the monopole 
defect in ISl E] , the direct physical interpretation of in the d = 2 case of the vortex defect is 
yet lacking. Both in the monopole and vortex cases parameter is involved into a condition 
for just one of the modes of the solution to the Dirac equation. 

When the transverse size of the vortex defect is shrinked to zero, the magnetic field 
strength takes form 

5(x) = 27r<l>5(x) , (2.7) 
and the vector potential can be chosen as 

y^(x) = -$r~^sin¥?, V"^( x) = $ r"^ cos . (2.8) 

Then Dirac Hamiltonian ()2.H) takes form 

H = -i-f'^Ydr - ir-SV(<9^ - ^e$) + 7°m , (2.9) 

where 

Y = 1'^ cos ip + •y'^ simp , 7*^ = — 7"*" sin + 7^ cos . (2-10) 

Using explicit form of 7 matrices (j2.3|) - (j2.4|) . one finds that the Dirac equation in the back- 
ground of a pointlike defect is reduced to following set of equations for the modes of Eq. ()2.5p : 

m dr + s{n - + s)r'^ \ ( ■^'^ \ = E ( ■^'^ ^ (2 11) 

-dr + s{n - e$)r-i -m J \ 9n J \ 9n J ' 

Partial Hamiltonians are essentially self-adjoint for all n, with the exception of n = uq, where 

no = le^ + ^-^s, (2.12) 

{u} is the integer part of quantity u (i.e., the largest integer which is less than or equal to 
u). Correspondingly, the modes with n ^ hq are regular at r = (i.e. at the location of 
the defect). The partial Hamiltonian at n = no requires a self-adjoint extension according 
to the Weyl-von Neumann theory of self-adjoint operators (see, e.g., Ref. j^), which upon 
implementation yields following condition for the corresponding mode ITHl IT^: 



cos 



+ ^) lim(|m|r)^/„o = -sgn(m) sin + lim(|m|r)^ ^g^o , (2.13) 
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where 

, . / 1, n > 
^g^^") = \ -1, u<0 

G is the self-adjoint extension parameter, and 

F = s{[e$]} + i - 1. , (2.14) 

{[■u]} = -u — [m] is the fractional part of quantity u, < {[m]} < 1; note here that Eq. ()2.1H|) 
implies that < F < 1, since in the case of F = - — -s both /^^ and obey the condition 

of regularity at r ^ 0. Note also that Eq. ()2.13p is periodic in B with period 2tt. 

So far solutions corresponding to the continuous spectrum, \E\ > \m\, are concerned, 
that obey the "orthonormality" condition 

d^x{E,n\^){^\E',n') = ^S^^£Xsnn' ■ (2.15) 

EE'\ 



Owing to Eq. ()2.1Hj) . an additional solution corresponding to the bound state with energy 
E = Ebs, \Ebs\ < \^\, appears at cos0 < 0, that obeys usual normalization condition 

d'^x{EBs,no\x){ x.\EBs,no) = 1. (2.16) 

Its energy is determined as a real root of algebraic equation 

(l+m-^EBsY~^ 

^ -A, (2.17) 



{l-m-^EBsY 
where 

T{u) is the Euler gamma function. The bound state energy is zero, Ebs = 0, at A = — 1; 
otherwise, we get 

sgn(FB5) = ^sgn(m)[sgn(l + A~^) - sgn(l + A)\ . (2.19) 

At cos© > (A > 0) the right hand side of Eq. ()2.19|) turns to zero, which corresponds to 
the absence of bound state in this case. 

Ill Resolvent and spectral density 

The kernel of the resolvent (the Green's function) of the Dirac Hamiltonian in the coordinate 
representation is defined as 

G-(r, ^; r', ^') = (r, - ^)->', ^'), (3.1) 
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where is a complex parameter with dimension of energy. The expansion of Eq. (j3.ip in 
modes takes form 



1 

2^ 



(3.2) 



In the case of H given by Eq. ()2.9p radial components of G^{r, ip; r', ip') ()3.2|) satisfy equations 
(compare with Eq. ()2.1H) ): 



-(uj + m) J V bnir;r') c„(r;r') 

dr' + s{n — + s)r'^^ \ f an{r;r') bn{r;r') 



— {uj — m) dr + s{n ~ + s)r an{r;r') dn{r;r') 

—dr + s{n — e^)r~^ 

— (CU — m) K^r' 

-dr' + s{n — e^)r'~^ —{uJ + m) J \ dnir;r') Cn{r;r') 



Off-diagonal radial components are expressed through the diagonal ones: 

bn{r; r') = {uj + m)^^ [—dr + s(n — e$)r^^] a„(r; r') = 

= (u — m)~^ dr' + s(n — + s)r'~^ 



c„(r;r'), (3.4) 



dn{r] r') = {uj — m) ^ \dr + s{n — e$ + s)r ^] c„(r; r') = 

= (uj + m)^^ ~dr' + s{n — e^)r' ^ an{r;r'). (3.5) 

In Appendix A we determine the diagonal radial components which can be presented in the 
following way, 

type 1 (/ = s{n — no) > 0): 



a„(r; r') = "-^{uj + m) e{r - r')H[]}p{kr)Ji^F{.kr') + e{r' - r)Ji_F{kr)H\'Jp{kr') , (3.6) 



r{i) 



Cn(r; r') = "-^{uj - m) \e{r - r')Hil\^p{kr)Ji+,Mkr') + e{r' - r)Ji+^_F{kr)H\'^,^p{kr') 



r(i) 



type 2 (/' = -s{n - Uq) > 0): 



(3.7) 



«n(r; = y (^ + m) [0(r - r')Hl^lp{kr)Ji'+pikr') + e{r' - r)Ji.+F{kr)H\^lp{kr')'^ , (3.8) 



r(i) 



ITT 



Cn{r; r') = y (cu-m) e{r - r')H[}\^p{kr)Ji'^i+F{kr') + e{r' - r)Ji' _i+F{kr)H\!'\^p{kr') 



r(l) 



(3.9) 
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type 3 (n = no): 

(^no{r;r') = ^ + ^ — j^i^^ _ r')H^_}l{kr)[smu^J^p{kr') + cos u^Jp{kr')] + 

2 sm + cos u^e^^^ I 

+ e{r' - r)[smu^J_Fikr) + cos u^JFikr)]H^^l{kr')^ , (3.10) 

Cno{r;r') = — — — \e{r - r')H[^}p{kr)[smiy^Ji_F{kr') - cosiy^J^i+F{kr')] + 

2 sm + cos z/^^e*-^ I 

+ e{r' -r)[smu^Ji^F{kr) - cosz/^J_i+F(A;r)]if(L^^(A;r')} . (3.11) 

Here k = \/ lJ^ — m^, Q{u) ~ ^ sgn(M)], J\{u) is the Bessel function of order A, Il^l\u) 
is the first-kind Hankel function of order A, and 

tanz/j^ = sgn(m) |m|^~^'^y4 , (3-12) 

u; + m 

where A is given by Eq. ()2.18|) . Note that the type 1 and type 2 components are regular at 
r = (or r' = 0), whereas the type 3 components are irregular at r = (or r' = 0), satisfying 
conditions (compare with Eq. (j2.13|) ): 

cos {^1^ + lim(|m|r)^a„(,(r; r') = — sgn(m) sin {^~^ + lim(|m|r)^~'^6„y(r; r'), 

(3.13) 

cos {^1^ + lim(|m|r)^(ino(r; r') = — sgn(m) sin {^1^ + lini(|m|r)"'^~'^Cno('"; ?"')? 

(3.14) 

and the ones at r' 0, which are obtained from Eqs. ()3.13|) and ()3.14|) by interchange 
Taking r' > r for definiteness, we get relations 

27r 

/ dyD tr G"^(r, r', y?) = ^ [a„(r; r') + c„(r; r')] = 

Q nGZ 

= ^ [(cu + m)Ii^F{-ikr)Ki^F{-ikr') + (u; - m)Ii+i^F{-ikr)Ki+i^F{-ikr')] + 



+ ^ [(u; + m)Iii+Fi-ikr)Kii+Fi-ikr') + {cu - m)Ii/^i+Fi~ikr)Ki/^i+F{-ikr')] + 



i'>i 



+ (a; + m)lF(-ikr)KF(-ikr') + (a; - m)Ii^F{-'ikr)Ki_F{-'ikr') H 2sin(F7r) ^ 

TT (tanz/^ + e*^'^) 



X 



[{u + m)ta.nu^KF{-ikr)KF{-ikr') + (tu - m)e^^^iri_ir(-zA;r)iri_^(-iA;r')] , (3.15) 
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and 



27r 



dip tr G"^(r, ip; r', ip) 



= ^ [an{r;r') + c„(r;r')] 



e<I>=0 



2lj In{-ikr)Kn{-ikr') , (3.16) 



where /a(w) is the modified Bessel function of order A, and 

TT 



2 sin(A7r^ 



is the Macdonald function of order A; note that the last equahties in Eqs. ()3.15|l and 
are obtained under condition Irak > 0. Using relations (see, e.g., Ref.jSTj) 



1 f dy 



2J y 



exp 



2y 2rr 



iy 
ey I due--h{u)~y{h{y) + h+. 



, Re\> -1, 



«6 
1>1 



we perform summation in Eqs. ()3.15|) and ()3.16p and get in the case of Irak > \Rek\: 

2 sin(F7r) 



2ir 

J dip tr 



G'-(r,<^;r',<^)-G-(r,<^;r',<^) 



e$=0. 



7r(tani^<^ + e'^"") 



[(cu + m) tan u^K F{i^r)K p{Kr') + {uj — m^je^^"" Ki_p(Kr)Ki_F{Kr' 

oo 

sin(F7r) f \ H?rr' (r — r')^ 



Lj dy exp 



21/ 



2rr' 



27rF(l - F) 



oo ^ 

(iMe-"[(l-F)ir^(n) + Firi_^(M)]-e-^(2F-l)[ir,.(l/)-iri_p(y)] I, (3.17) 
y ) 

where k = —ik. Taking the limit r' — r in Eq. ()3.17|) and integrating it over the radial 
variable, we get the renormalized (finite) trace of the resolvent operator 



oo 2lT 



Tr (H — uj) ^ = J drr J diftr 







G'^{r,!p;r,!p) - G^{r,ip;r,ip) 
2F-1 



e<S>=0. 
1 



note that the last result can be continued analytically to the whole complex cu-plane. Note 
also that Eq. (j3.18p can be rewritten in an equivalent form: 



Tr {H -uo)'^ 
1 



^2 _ ^2 



2F -1 



g-i(l-F)7r u^-l 



[l-Ff 



+ m 



1 



+ 1-F 



(3.19) 



Taking the imaginary part of Eq. ()3.18p or Eq. ()3.19|) at u = E + iO,we get spectral density 
t{E), see Eq.dEH). 



IV Thermal average and fluctuation of the charge 

Taking into account Eq. (jl.9j) . one can get the following contour integral representation for 
induced charge (jl.8p and its quadratic fluctuation (jl.l(J|) : 



and 



c 



c 



(4.1) 



(4.2) 



where C is the contour (— oo + iO, +oo + iO) and (+oo — iO, — oo — iO) in the complex cu-plane. 
Substituting Ea. dTTHIl for Recu > and Eq.(jSIIll) for Recu <0 into Eqs.(jiIIl) and (jOJ, we 
obtain 



QiT) 



sgn(m) 1^ [sgn (l + A'^) - sgn(l + A)] tanh (^^P\Ebs\^ + 

oo 

2sin(F7r) f du , /^l ^, , / 

^ ^ ' =tanh (-P\m\Vu + lj x 



vr J Uy/u + 


,F n IP\ A-l„,l~F 



and 



A2 



FAu^ - (1 - F)A-^u^~^ - u cos(F7r) + (F - i) u {Au^ + A-^m^-^) 
[Am^ - + 2 cos(F7r)]^ + A{u + 1) sin2(F7r) 



^ |^[1 - sgn(A)] sech^ (\f3\EBs\) - i^(l - F) sech^ Q/5|m| ) + 



(4.3) 



+ 



2 sin(F'7r) /" c/m 



TT 



y ^ sech^ Q/5|m|VM + 1^ 



FAm^ + (1 - F)^-!^!-^ - (2F - l)ucos(F 



[Am^ - A-iyi-^ + 2 cos(F7r)]' + 4(m + 1) sin2(F7r) 



(4.4) 
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where A is given by Eq. (j2.18|) . 

In the cases oi A = and A~^ = expressions for the charge and its fluctuation simphfy: 



Q(T) = - i ± M tanh Q/3mj , 6 = ±s^ (mod27r), 



and 



A2 



, F - - ± sech^ (-P\m 
4 V 2 2/ 12^^' 



TT 



e = ±s- (mod27r); 



note that Eg. ()4.5|1 at O = s— (mod27r) was obtained in Ref.jZj. 

In the hmit T ^ (/3 ^ oo) the charge tends to finite value (see Ref.|19j): 



(4.5) 
(4.6) 



g(o) 



-sgn(m)(l -F), 

g 

- - sgn(m)F, 

g 

- -sgn(m)(l + F), 

- - sgn(m)F, 

g 

-sgn(m)(l - F), 

g 

-sgn(m)(2 - F), 



yl 



-1 < A < oo 
-1 = -1, = 
-oo < A< -\ 

-1 < A'^ < oo 
A = -1, A = 
-oo < < -1 



> , < F < 



< F < 1 



(4.7) 



Q{0) 



0, 



— s sgn(m) arctan ^tan -^^ , 7^ vr (mod 27r) 

= TT (mod27r) 



(4.8) 



whereas the fluctuation tends exponentially to zero for almost all values of B with the 
exception of one corresponding to the zero bound state energy, Ebs = {A = —1): 







'^0(0) 



2 

-1 



4 ' 



(4.9) 



In the high-temperature limit the charge tends to zero: 



e , ,sin(F7r)r(l-F) / 9 tt 
- sgn(m) — : — ^i^^ tan | s— + — 



vr r(l + F) 



nfi- \ J e sm 
Q{T ^ 00 = <^ - g sm e, 
8 Kb-I 



\m\ 



1-2F 



e , ,sin Fvr T F / 6 7r\ / m 

-- sgn(m) — ^ , , cot s h - — - 

2 ^ ^ ' 71 r(2-F) V 2 47 \kBT 



2F-1 



o<F<l 

2 

- < F < 1 
2 



(4.10) 
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whereas the fluctuation tends to flnite value, see Appendix B: 



n 



-{l-F)\ e^s-{mod2n) 



2 

—F 
4 ' 



e = s- (mod27r) 



TT 



e 7^ -s- (mod27r) 
2 



TT 



— {l-FY, e = -s-(mod27r) 



At half-integer values of e$ one has 



= tan I s— + - 



and the charge and its fluctuation take form 



— |[1 — sgn(cos©)] tanh ^(3msmQ^ + 



+ 



sin 29 



dv 



> , < F < 



< F < 1 



tanh ( -/9mA/u 



27r y ^/^(v - 1) V - sin^ © 



(4.11) 



(4.12) 



(4.13) 



and 



I 



^ \ [1 — sgn(cos©)] sech^ [ i/?|msin©|j — scch^ f ) + 



+ 



cos Of dv 

TT J y/v-1 



sech^ ^^/9|m|v^^ 
V - sin^ © 



(4.14) 



An alternative representation for the charge and its fluctuation is obtained by deforming 
contour C to encircle poles of the tanh ^^/^"^^ sech^ ^^'^'^^ functions, which occur 

cUn — (2n + J •> Appendix C: 



Q{T) = -esgn(m) <! " ^ tanh (^^ | + 



e V- 2(2F- l)(2n + 1)2^2 + A[l + (2n+ 1)2^2]-^ -A-i[l + (2n + l)2e 



TT ^ [l + (2n+l)2e2]{/l[l + (271 + 1)2^2]^ + 2 + ^-i[l + (2n+l)2^2ii-Fi f 

ng Z 
n>0 



(4.15) 
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and 



AJ(T) = ^[1-2F(1-F)]sech2 (^^) + 



n>0 

X |(2F - l)[(2n + 1)2^2 - + {2n + l)2eY - ^^'[1 + (2n + 1)'^']'"''}+ 

+ 2{1 - [3 - 4F(1 - F)]{2n + l)^^^! _ 4(2^ ^ 1)2^2^ 
(2F - l){v4[l + {2n + 1)2^2]^ _ ^ + _ ^ + (2F - l)2(2n + 1)2^2 

^ A[l + (2n + 1)2^2]^ + 2 + + (2n + 1)2^2]!-^ 

(4.16) 

where ^ = 7r/(/3|m|). 

V Discussion 

In the present paper we consider an ideal gas of twodimensional relativistic massive electrons 
in the background of a static pointlike magnetic vortex. This system at thermal equilibrium 
is found to acquire electric charge: its average Q{T) is given by Eq. ()4.3|) . and its quadratic 
fluctuation ^^q(t) is given by Eq. ()4.4|) . The most general boundary conditions (parametrized 
by the self-adjoint extension parameter B) at the location of the vortex are employed, and 
arbitrary values of the vortex flux $ are permitted; our results are periodic in G with period 
27r at fixed $ and periodic in $ with period at fixed B (e is the electron charge). 
Note that Eqs. ()4.3|l and ()4.4|l can be regarded as the Sommerfeld- Watson transforms of the 
infinite sum representation, Eqs. ()4.15j) and ()4.16j) . Note also that the charge is odd and its 
fiuctuation is even under transition to the inequivalent representation of the Clifford algebra 
(s —s or m ^ — m). 

Eq. ()4.3|) can rewritten in the form 

g(T)=Q(0) + Q(T), (5.1) 

where (5(0) is given by Eqs. (pT7jl - (prK|l [TTI]. and 

~ e /sgn(l + A-i) -sgn(l + .4) , 2F - 1 

Q[T) = -sgn(m) ^ v / , 



2 " ^ ^ i exp (PIEbsD + 1 exp(/5|m|) + 1 



fj 1,2 A 



+ — — / dw sech. ( -/3|m|w ) arctan 



A{w^ - If - A~\w^ - ly-^ + 2 cos(F 



1 



2w sin(F7r) 



(5.2) 



recall that F is related to the fractional part of e$ by Eq. ()2.14|) . A is related to O by 
Eq. ()2.18j ). and bound state energy Ess is determined implicitly by Eq. ()2.17p . 
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Our result should be compared with the result of Refs.j2Sl 1211 123 




) 



(5.3) 



where $ is the flux of a magnetic fleld with an extensive support, and it is implied that the 
region of the support is not excluded. Thus, result ()5.3|1 describes the direct effect of the 
field strength, whereas our result describes the indirect, through the vector potential, effect 
of the field strength from the excluded region. In contrast to Eq. ()5.3p . our expressions for 
Q(T) and ^q{t) periodic in the value of the fiux, vanishing at integer values of e$, and 
this can be regarded as a manifestation of the Bohm-Aharonov effect in quantum field 
theory at nonzero temperature. 

The nonvanishing of the charge quadratic fiuctuation signifies that the charge of the 
system is not a sharp quantum observable and has to be understood as a thermal expectation 
value only. In the high-temperature limit the average charge tends to zero (j4.1Uj) and the 
fiuctuation tends to finite value ()4.11|) . In the zero-temperature limit quantities ^q(^t) 
Q{T) tend exponentially to zero and the charge becomes a sharp quantum observable with 
finite value Q{0) ()4.7|) - ()4.8p . However, the last statement is true for almost all values of 
with the exception of one corresponding to the zero bound state energy, Ebs = (A = — 1), 
since in this case the zero-temperature fiuctuation is nonzero, see Eq. ()4.9|1 . 

At half-integer values of e$ the average charge takes form of Eq. ()4.13|l which coincides 
(after substituting s for 2eg, where g is the magnetic monopole charge, 2eg = n is the 
Dirac quantization condition) with the expression for the thermally induced charge in the 
monopole background in threedimensional space [Zl El E] • It should be emphasized that at 
non-half-integer values of e$ the behaviour of the charge as a function of differs drastically 
from the one at half-integer e$. 

To see this explicitly, we plot the charge and its fluctuation as functions of for several 
values of the vortex flux and temperature on Figs|T]|ni F = 0.1, 0.3, 0.5, 0.7, 0.9. Here 
values {kBT/\m\) = 5~\ 1, 5 correspond to two dashed (with longer and shorter dashes) 
and one dotted lines, and values T = and T = oo correspond to solid lines; the latter 
cannot lead to confusion, since, as it has been already noted, the charge at T = oo equals to 
zero everywhere, while the fluctuation at T = equals to zero almost everywhere with the 
exception of one point {A = —1). Two dashed lines coincide practically in the utmost left 
parts of Figs, la, 2a and in the utmost right parts of Figs. 4a, 5a. The qualitative difference 
between the F = 1/2 and F ^ 1/2 cases is most evident at zero temperature and is persisting 
with the increase of temperature, notwithstanding the dying of the charge on the whole at 
high temperature. 

In the 7^ 2 charge at zero temperature is given by a step function with two 

jumps. As temperature increases, the jump corresponding to the zero bound state energy 
[A = —1) is smoothed out, while another jump is persisting. The charge at A = —1 is not a 
sharp quantum observable even at zero temperature, which is explicated by the nonvanishing 
of the fluctuation in this case. As temperature departs from zero, the fluctuation develops 
a maximum at A = — 1 and a minimum close to the position of the persisting jump of the 
charge, but out of the region where bound state exists. With the increase of temperature the 
maximum is widening and disappearing, while the minimum is narrowing with its position 
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approaching the position of the charge jump and its width tending to zero in the high- 
temperature hmit. 

In the ^ ~ 2 charge at zero temperature is hnear in with one jump at A = —1 

(0 = s7r(mod27r)) where the charge is not a sharp quantum observable. As temperature 
increases, this jump is smoothed out. Appropriately, the fluctuation is symmetric with 
respect to the position of this jump, and a maximum of the fluctuation is smoothed out with 
the increase of temperature. 

In conclusion we note that the system considered can acquire, in addition to the charge, 
also other quantum numbers. In the case of zero temperature this issue is comprehensively 
elucidated in Refs.jini IHB]; and an appropriate generalization to the case of nonzero 
temperature will be studied elsewhere. 
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Appendix A 

The diagonal elements of G"^(r, r', </?') (j3.2j) satisfy second-order equations: 

[-r^'^drvdr + r'^{n - e$)^ - + m^] an{r\ r') = {u + m) ^^^^^\ (A.l) 

[-r~^drrdr + r~^{n - e$ + sf - cu^ + m^] c„(r; r') = (to - m)-—j=^. (A.2) 
The general solution to, say, Eq. ()A.l|) has the form 

a„(r; r') = ^(o; + m) [e{r - r') [Hfl_^^^{kr)Hfl_^^^{kr')- 

-<Le.)(^0<le.)(^^')] + <Le.)(^Opi"nr') + ^Le.) (^OptH^') } , (A.3) 

where 

Hx\^) = -4v-^ [e"'"^'^A(«) - J-x{u)] and H^^\u) = [J^x{u) - e^'^Uu)] 

are the first- and second-kind Hankel functions of order A. Without a loss of generality one 
can choose a physical sheet for square root k = ^uS^ — vn? as < Arg k < ti [Imk > 0). 
Then we impose the condition that solution ()A.3|) behaves asymptotically (at r — > oo) as an 
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outgoing wave ^ ^ J , and this yields: pn {r') = H^^^_^^^{kr') . Thus we get 



r — r 



(2) 



+^(^' - r)H%_^^pr)H%_^^^{kr')\ + <l,^)(A:r)p(r)(r')} . (A.4) 
In a similar way we get for the solution to Eq. ()A.2|l 



+^(-' - + ^f;L*)+i(A:r)pi'=)(r')} . (A.5) 



Quantities p"n\r') and pn'{r') are determined by the condition at r — > 0. As it has been 
discussed in Section 11, the condition of regularity at r — > is imposed in the case of n 7^ ng, 
and this yields: pn\r') = Hf^^_^^-^{kr') and pn {r') = Hlll_^^)+i{kr') {n ^ hq). Thus we 
get 



a„(r;r') 



~2 



+ m) e{r - r')H\l^_^^^{kr)J\^_,^\{kr') + 



+e{r' -r)J\^_,^\{kr)H\^_^^^{kr') , n ^ n^, (A.6) 



and 



Cn{r; r') = —{u - m) ^9{r - r')/ff„'i,$+,|(A;r) J|„_e$+«|(A;r') + 

'I 



n ^ no, (A.7) 



which gives the type 1 and the type 2 solutions (|3.6|) - (j3.9|) . 

In the case of = no the solutions to Eqs. ()A.l|) and ()A.2|) are not regular at r — 0, 
but their irregular behaviour has to be matched with the one of the 6„o ^"^^ '^no components 
correspondingly, owing to conditions ()3.13p - ()3.14p . Using Eqs. ()3.4|) and ()3.5|) . we get 



e{r - r')H[%{kr)H^yp{kr') + 

+eir' - r)H^i'}p{kr)H^yj,{kr')j + H[%{kr)pl^^{r')^ , (A.8) 



d„^(r; r') = { d{r - r')H%{kr)Hflp{kr') + 

+6{r' - r)H^%{kr)Hflp{kr')\^ + H%{kr)p'^l{r')^ . (A.9) 

Substituting the pair of Eq. ()A.4|) at n = no and Eq. ()A.8|) into Eq. ()3.13p and the pair of 
Eq. ()A.5|) at n = rio and Eq. ()A.9|) into Eqs. ()3.14|) . we determine Pn}{r') and pnl{r'), and 
obtain the type 3 solution ()3.10p - ()3.1H) . 
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In the absence of the vortex defect radial components for all n are regular at r — > 0: 



+e{r' - r)Jn{kr)H^^\kr')\ , (A.IO) 



= '-^{uo - m) [e{r - r')H^2\kr)Jn{kr') + 



+d{r' - r)Jn{kr)Hi^\kr')] . (A.ll) 



Appendix B 

In the high-temperature limit Eq. ()4.4|) takes form 



hm 

T^oo 



.2 n 



2 sin(F7r) f du FAu^ + (1 - F)A-^u^~^ - (2F - l)u cos(F 



TT 



[Am^ - A-^u^-^ + 2 cos(F7r)]^ + 4(m + 1) sin2(F7r) 



(B.l) 



Using relation 

d {Au^ + A'W-^) sin(F7r) 

- A-^u^-P) cos(F7r) + 2 ~ 
_ 2 sin(F7r) FAm^ + (1 - F)A-^u^-^ - (2F - l)u cos{Ftt) 
u [Au^ - A-^u^-P + 2cos{F7r)f + 4:{u+l)sin^{F7r) 

we get in the case of — oo < A < 0: 



(B.2) 



2 sin(F7r) f du FAu^ + (1 - F)^-^^^-^ - (2F - l)u cos(F7r 



TT 







M [v4m^ - A-^u^~^ + 2 cos(F7r)]^ + 4(u + 1) sin2(F 



1 {Au^ + A-i^i-^) sin(F7r) 

— cLrCtcin — ; = —— — ; 

TT {Au^ - A~^u^-P) cos(F7r) + 2 



u=0 



1 (AM^ + A-V-^)sin([(l-F)7r] 

^ '^^''^'^^ {-Au'' + A-ini-^) cos[(l - F)u] + 2 



-F 



u=0 



< F < 



-1 + F, 2<^<1 



(B.3) 



In the case of < A < oo one should note that the integration extends over both the 
principal and neighboring sheets of the Arctan function (here Arctan z = arctan2; + n7r, and 
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tan(Arctan z) = z): 

OD 

2 sin(F7r) f du FAu^ + (1 - F)A-^u^-^ - (2F - l)u cos(F7r) 



TT 







u [AuP - A-^u^-P + 2 cos(F7r)]^ + 4(u + 1) sin2(F7r) 



1 (Au^ + A-^mI-^) smfFvr) 

1 H — arctan — - — = — — - — ; 

TT {Au^ - A-H^-P) cos(F7r) + 2 

1 + A-V-^)sm([(l -F)7r] 

+ - arctan ^ A-i^i-^) cos[(l - F)7r] + 2 



M=0 



M=0 



< F < - 
2 



< F < 1 



(B.4) 



In the cases of A = and A ^ = we use Eq. ()4.6|) . and in the case oi F = - we use 
Eq. ()4.14p . Thus, we get Eq. ()4.11|) as the high-temperature hmit of the fluctuation. 



Appendix C 

If a function of complex variable uj has a pole of Z— th order at a; = Un-, then the integral 
over a contour encircling this pole is given by expression 

du fiuj) = 2m\im^ (7^^ [/(^)(^ " ^")'] • 

By deforming contour C in Eqs. fj4.H) and ()4.2|) to encircle poles of first and second orders, 
we get 

QiT) = -^Y.Tr{H-u;r.)-\ (C.l) 

and 

^lm = -j,J2^r{H-u;^)-\ (C.2) 



nS 1 



I'jl _____ 

where a;„ = (2n + 1)-^- Using Eq. ()3.18|) . we get 

P 



Tr {H - uo)-^ = — <j F{Fuj^ + 2ujm + Fm"^ 



^^ {2F-l)uj' + 2um + m' [(2F - 1 V + m]^ 

e-^^- tan + 1 [g-*^- tan z/^ + 1]^ ' ' 



Consequently, we obtain 



1 + (2F - l)i(2n + 1)^ sgn(m) 



1 + A[l + (2n + + ^(^2n + 1)^ sgn(m)]-i 



(C.4) 
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and 



l-(2F-l)(2n + l)^e' + 2z(2n + l)esgn(m) 
1 + A[1 + {2n + l)2^2]F[i ^ ^(2n + l)^sgn(m)]-i 

[l + (2F-l)2(2n + l)esgn(m)]2 ) 
{l + A[l + (2n + l)2e2]^[l + 2(2n + l)esgn(m)]-i}2j ' ^ 

where ^ = n/{P\m\). Summing over n, we get Eqs. (j4.15|) and (|4.16|) . 
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Figure 1: F = 0.1 
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Figure 2: F = 0.3 
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Figure 3: F = 0.5 
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Figure 4: F = 0.7 
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